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Abstract 

Pleurotus eryngii has recently become a major cultivated mushroom; It uses tetrapolar heterothallism as a part of its 
reproductive process. Sexual development progresses only when the A and B mating types are compatible. Such mating 
incompatibility occasionally limits the efficiency of breeding programs in which crossing within loci-shared strains or 
backcrossing strategies are employed. Therefore, understanding the mating system in edible mushroom fungi will help 
provide a short cut in the development of new strains. We isolated and identified pheromone and receptor genes in the S3 
locus of P. eryngii and performed a functional analysis of the genes in the mating process by transformation. A genomic 
DNA library was constructed to map the entire mating-type locus. The B3 locus was found to contain four pheromone 
precursor genes and four receptor genes. Remarkably, receptor PESTE3.3.1 has just 34 amino acid residues in its C-terminal 
cytoplasmic region; therefore, it seems likely to be a receptor-like gene. Real-time quantitative RT-PCR (real-time qRT-PCR) 
revealed that most pheromone and receptor genes showed significantly higher expression in monokaryotic cells than 
dikaryotic cells. The pheromone genes PEplibS.l and PEplibS.S and the receptor gene PESTE3.3.1 were transformed into P5 
{A3B4). The transformants were mated with a tester strain (A4B4), and the progeny showed clamp connections and a normal 
fruiting body, which indicates the proposed role of these genes in mating and fruiting processes. This result also confirms 
that PESTE3.3.1 is a receptor gene. In this study, we identified pheromone and receptor genes in the B3 locus of P. eryngii 
and found that some of those genes appear to play a role in the mating and fruiting processes. These results might help 
elucidate the mechanism of fruiting differentiation and improve breeding efficiency. 
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Introduction 

Pleurotus eryngii is a white rotter of wood and the only species 
of its genus that grows on the living stems and roots of umbellifers 
[1]. Pleurotus sp. mushrooms have rapidly emerged as one of the 
most popular cultivated edible mushrooms [2]. In particular, P. 
eryngii, the king oyster mushroom, has recently become a major 
cultivated mushroom in Korea, China, and Japan because of its 
good taste, flavor, and nutritional factors. As the consumption of 
P. eryngii increases, new strains need to be bred with improved 
traits. However, mating incompatibility occasionally limits the 
efficiency of breeding programs in which crossing within loci- 
sharing strains or backcrossing strategies are employed. When a 
specific strain is required to have a single useful trait such as 
disease resistance, backcrossing is considered to be the best 
method for introducing the trait into the recurrent parent. In a 
backcross, however, only 25% of all possible combinations are 



compatible in tetrapolar mushroom fvmgi. Therefore, a breeder 
cannot evaluate the traits of incompatible hybrids. Although the 
mating-type systems of several edible mushrooms have been 
examined [3,4], most current knowledge has been derived from 
studies of model mushrooms such as Schizophyllum commune and 
Coprinopsis cinerea [5,6]. Thus, the mechanisms controlling the 
mating system in P. eryngii require further examination to 
improve the efficiency of P. eryngii breeding programs and to 
address its basic biology. 

P. eryngii uses tetrapolar heterothallism as a part of its 
reproductive process. In a tetrapolar system, genes in the A locus 
encode a homeodomain transcription factor, whereas the B locus 
encodes both peptide pheromones and receptors (for review, see 
[7]). The structure of the B locus differs in different species. For 
example, in Ustilago maydis, a single locus contains one receptor 
and one pheromone, whereas S. commune possesses two 
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functionally independent sub-loci containing one receptor and 
several pheromones [6]. In C. cinerea, a common structure is three 
sub-groups, with one receptor and multiple pheromones in each 
sub-group [5,8,9]. However, recent genome sequence reports have 
shown that the matB locus is much more complex than previously 
realized [4,10,1 1]. Some receptors referred to as receptor-like were 
found at loci distinct from the known mating loci, and until now, 
whether the receptor-like gene(s) had a role in the mating process 
was unclear. 

Previously, sixteen alleles at locus A and 15 alleles at locus B 
were identified from the 12 strains of P. eryyigii, and the B3- 
specific SCAR primers discriminating the B3 locus were 
developed based on the RAPD (random amplified polymorphic 
DNA)-derived sequence, which was specifically amplified from B3 
locus-harboring strains [12]. In this study, we isolated the mating- 
type B3 locus of P. eryngii using a PCR-based library screening 
method with a specific marker. A physical gene map of the region 
controUing mating type with both flanking regions is presented 
here. The expression patterns of the pheromone and receptor 
genes were analyzed in both monokaryotic and dikaryotic mycelia. 
Functional analyses of the plun-omonc and receptor genes were 
performed hy AgrohacleriiimArcmcd transformation. To the best of 
our knowledge, this study is the first functional analysis of 
pheromone and receptor genes in edible mushrooms. 

Materials and Methods 

Fungal strains and growth conditions 

Two monokaryon strains, P5 and P6 (Table 1), were obtained 
from the dikaryotic mushroom KNR2312 using de-dikaryotization 
through protoplast homogenization [12] and has been deposited at 
Korean Culture Center of Microorganisms as KCCM90103 (P5) 
and KCCM90104 (P6). Cultured mycelia were obtained by 
inoculating 3—4 pieces of 1 x 1 cm^ MCM (mushroom complete 
medium: 0.2% peptone, 0.2"'<) yeast extract, 2.0% glucose, 0.05% 
MgS04-7H20, 0.05% K2HPO4, and 0.046% KH2PO4) agar with 
mycelia in MCM broth, followed by incubation for 14 days at 
25°C and 120 rpm. After incubation, the myceUa were filtered 
with sterile Whatman filter paper No. 41 and lyophUized before 
grinding with a sterile mortar. For long-term storage, strains were 
cultivated on slant medium of MCM, filled with sterile mineral oil 
(Sigma, St. Louis, MO, USA), and placed at 4°C until further use. 

Fosmid library construction and screening 

High-molecular-weight P6 DNA was purified using the 
Neurospora SDS buffer method as described previously [13], with 
minor modifications. First, 1 g of lyophilized mycelium powder 
was added slowly to 10 mL oi Neurospora SDS buffer (0.15 M 
NaCl, 0.1 M EDTA, and 2% SDS, pH 9.5) in a 125-mL flask. 
Second, 0.5 mL of 2 mg/mL proteinase K (Fermentas, Waltham, 
MA, USA) solution was added, followed by incubation at 37°C for 
up to 24 h with gentle agitation; 10 mL of sterile water was added 
to the slurry, and the cellular debris was removed by centrifugation 
at 1 5,344 X g. The supernatant was extracted 3-5 times with 
phenol saturated with Tris-HCl (pH 8.0) and once with water- 
saturated chloroform. DNA was precipitated by adding 2.5 
volumes of absolute ethanol and removed with a sterile glass rod 
into a 1.5-mL centrifuge tube. After briefly air-drying, 1 mL of 
sterile TE buffer (10 mM Tris-EDTA, pH 8.0) containing 30 |xg/ 
mL RNase A (Qiagen, Hilden, Germany) was added. After 
incubation at 37°C for several hours with occasional mixing by 
inversion, the genomic DNA (gDNA) solution was stored at 5°C 
until use. The gDNA was used to construct a library using the 
CopyControl Fosmid library production kit (Epicentre Biotech- 



nologies, Madison, WI, USA) according to the manufacturer's 
instructions. 

Screening of clones harboring the mating-type region from the 
fosmid library was conducted using direct PCR; 25 clones were 
pooled in LB+chloramphenicol (12.5 |J,g/mL) using sterile tooth- 
picks, followed by incubation for 24 h at 37°C and 120 rpm. 
Plasmid DNA was extracted with the GeneAll EXPREP plasmid 
SV MINI Kit (GeneAll Biotechnology, Seoul, Korea) according to 
the manufacturer's protocol. The 53-specific SCAR primer set 
(Table SI) was used to specifically amplify strains harboring the B3 
locus. The thermal cycling reaction was conducted in 20-|j,L 
volumes using GoTaq (Promega, Madison, WI, USA) in a Dyad 
thermal cycler (Bio-Rad, Hercules, CA, USA) with the following 
parameters: initial denaturation at 94°C for 4 min, followed by 35 
cycles of denaturation at 94°C for 1 min, annealing at 65°C for 
1 min, and extension at 72°C for 1 min 30 s. Two clones carrying 
the B3 locus were isolated by further PCR with individual clones 
displaying the SCAR marker band. 

Fosmid sequencing and data analysis 

To isolate fosmid DNA, a QIAprep 96 Turbo Miniprep kit 
(Qiagen) was used. Two fosmids containing mating type-specific 
regions were sequenced by the shotgun method. Purified DNA was 
sheared into random fragments. After end repair and size 
fractionation, a DNA fragment was hgated into pUCl 18/i/mcII 
(Takara, Shiga, Japan) and then sequenced using the Big-Dye 
cycle sequencing kit and an ABI3700 or 3730x1 DNA analyzer 
(Applied Biosystems, Foster City, CA, USA). Base calling, contig 
assembly, and editing were conducted using the phred/phrap/ 
consed package (University of Washington, USA). 

Genes in the mating type-specific region were identified using 
DNAMAN (Lynnon Biosoft, Quebec, Canada) and FgeneSH 
(Softberry, http://linuxl.softberry.com/berry.phtml) with a C. 
cinerea model; the BLASTX function of the NCBI database was 
used to identify open reading frames (ORFs) and genes. ORFs 
with small peptide sequences (10-100 amino acids) were manually 
curated to identify the CaaX motif in the C-terminus and were 
compared with known sequences in the NCBI database using 
BLAST to identify small genes such as those encoding phero- 
mones. Additionally, the characterization and identification of 
putative mating-type genes were manually curated using the NCBI 
BLASTX and Protein BLAST databases (USA). Transmembrane 
motif analysis of receptor genes from P. eryngii and other receptor 
sequences retrieved from GenBank were predicted using PHO- 
BIOUS (http://phobius.sbc.su.se/), HMMTOP (http://www. 
enzim.hu/hmmtop/), and TMHMM (http://www.cbs.dtu.dk/ 
services/TMHMM/) software. Phylogenetic analysis was per- 
formed according to a previous report [4] . The receptor sequences 
of S. commune, P. djamor, C. cinerea, and U. maydis, available in 
the NCBI GenBank database, were used with previously reported 
accession numbers [14]. The receptor sequence oiLaccaria bicolor 
was obtained from a genome site at the Joint Genome Institute 
JGI: http://genome.jgi-psf org/Lacbi2/Lacbi2.home.htnil) ac- 
cording to an accession number listed in a previous report [15]. 
The U. maydis PRA2 was used as an out-group. Protein sequences 
were aligned using ClustalW 1.64, and maximum likelihood 
analyses were conducted to estimate phylogenetic relationships 
using RAxML 7.04. Branch support was inferred from 1000 
repetitions of nonparametric bootstrapping. 

RNA extraction and cDNA synthesis 

Total RNA from KNR2312 (dikar>'on) and P6 (monokaryon) of 
P. eryngii were extracted from lyophiUzed mycelia using modified 
protocols for TRIzol (Invitrogen, Carlsbad, CA, USA) and the 
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Table 1. Pleurotus eryngii strains used in this study and their mating interactions. 





Species 


Description 


Note (reference)^ 


KNR2312 


dikaryon, P5 {A4B3)+P6 iA384) 




P5 


monokaryon, A3B4 


protocione [12] 


P6 


monokaryon, A4B3 


protoclone [12] 


KNR2312-4 


monokaryon, A3B3 


meiotic spore [12] 


KNR2312-6 


monokaryon, A3B4 


meiotic spore [12] 


KNR23 12-29 


monokaryon, A4B3 


meiotic spore [12] 


KNR2312-2 


monokaryon, A4B4 


meiotic spore [12] 


P5phb31 


monokaryon, A3B4::PEphb3.1 


transgenic, this study 


P5phb33 


monokaryon, A3B4::PEphb3.3 


transgenic, this study 


P5STE331 


monokaryon, A3B4::PESTE3.3.1 


transgenic, this study 


P5phb31 XKNR2312-2 


dikaryon, positive control 


this study 


P5phb33xKNR2312-2 


dikaryon transformation 


this study 


P5STE331 XKNR2312-2 


dikaryon transformation 


this study 


P5xKNR2312-2 


negative control 


this study 



^Protoclones were prepared from the homogenization of protoplast [12]; meiotic spores were derived from meiotic recombination. 
doi:l 0.1 371 /journal.pone.Ol 04693.1001 



RNeasy Mini Kit (Qiagen). After extracting total RNA, the 
RNase-Free DNase set (Qiagen) was used to eliminate gDNA 
contamination. Both sets of RNA were evaluated for quantity and 
quality using a NanoDrop Spectrophotometer ND-1000 (Nano- 
Drop Technologies, Wilmington, DE, USA). The remaining 
samples were flash-frozen in liquid nitrogen and stored at — 70°C 
until use. Then, the isolated RNA was transcribed into cDNA 
using a reverse transcription system with oligo-dT primer 
(Promega). The converted cDNA was used for both sequence 
verification and gene expression analysis. Four putative phero- 
mone precursor genes and 4 receptor genes were amplified from 
the P6 cDNA using specific primers (deposited with the cDNA 
sequences in GenBank). PGR products purified from the gel using 
the QIAquick PGR Purification Kit (Qiagen) were transformed 
into Escherichia coli plasmids using the TOPO cloning kit for 
sequencing (Invitrogen) and the pGEM-T Easy Vector (Promega). 
Plasmids in transformed E. coli were isolated using Hybrid-Q 
(GeneAll Biotechnology) and prepared for sequencing by the 
GenoTech Gompany (GenoTech Gorp., Daejeon, Korea). After 
qualifying the cDNA sequences, the four pheromone and four 
receptor gene sequences, along with their primers for amplifica- 
tion, were deposited in GenBank under the following accession 
numbers: KG879316 for PESTE3.3.1, KG879317 for 
PESTE3.3.2, KG879318 for PESTE3.3.3, KG879319 for 
PESTE3.3.4, KG879320 for PEphb3.1, KC879321 for 
PEphh3.2, KG879322 for PEphb3.3 and KG879323 for 
PEphh3.4. 

Analysis of gene expression by real-time qRT-PCR 

Real-time qRT-PGR (GFX96, Bio-Rad) was performed using 
KNR2312 (dikaryon) and P6 (monokaryon) to quantify the 
expression of the 4 receptor genes and the 4 pheromone genes. 
The primer pairs used for real-time qRT-PGR amplification of the 
internal fragments of these genes are listed in table S2. Synthesized 
cDNA (100 ng/ ul) was used as a template for real-time qRT-PGR. 
Real-time qRT-PGR was run as a duplex, in which one duplex 
partner was used as the internal standard gene (18S rRNA of P. 
eryngii), with the accession number FJ572254. The other duplex 
partner was one of the 4 pheromone genes or 4 receptor genes. 



The Rotor-Gene Q, (QIAGEN, Hilden, Germany) was used for 
the reaction as follows: an initial step at 95°G for 5 min, followed 
by 39 cycles at 95°G for 5 s and 62°G for 10 s. One microlitre of 
cDNA was amplified in a 25.0 [xl reaction using Rotor-Gene 
SYBR Green PGR kit master mix (Qiagen) with each primer at a 
final concentration of 500 nmol/1. The Gt value of each gene was 
normalized for differences in the amount of total cDNA in the 
reaction using the 18S rRNA of P. eryngii as an internal standard 
control. The normalized values for the pheromones and receptors 
were then calibrated to the value for the monokaryon (P6), which 
was assigned a value of 1, by the standard curve method using the 
Rotor-Gene software (QIAGEN). All assays were performed in 
triplicate. The statistical analysis of the gene expression was 
performed by one-way analysis of variance (ANOVA) and Tukey's 
test (alpha = 0.05) using the SAS program (SAS 9.1, SAS Institute 
Inc., Gary, NG, USA). 

Agrobacterium-mediated transfornnation 

Genomic fragments containing a 1-kb upstream fragment for 
the promoter and a 0.7-kb downstream fragment for the 
terminator of two pheromone genes (PEphh3.1 and PEphb3.3) 
and one receptor gene [PESTE3.3.1) were ligated into the Sal I 
site of the pBGgHg binary vector carrying the hygromycin B 
resistance gene {hph) [16]. Although this vector has the Agaricus 
bisporus gpdil promoter for expression, we adapted each gene's 
own promoter and terminator. Recombinant pPEphbS.l, 
pPEphb3.3 and pPESTE3.3.1 were introduced into Agrohacter- 
ium tumefaciem AGL-1 and transformed into P. eryngii P5 as 
described previously [17], with minor modifications. Agrobacter- 
ium induction medium (IM) culture and mycelium (harvested from 
colonies grown on MGM and homogenized for 30 s) culture 
(100 |J.L each) were co-cultivated, followed by treatment with 
sonicatioii for 3 min and vacuum for 30 min, and then inoculated 
overnight at 25°G in darkness. 

Putative transgenic lines showing hygromycin resistance were 
confirmed by PGR using specific primers (Table SI) as described 
previously [17]. The SSR45 primer set showing co-dominant 
polymorphism on KNR2312, P5 and P6 [18] was used to confirm 
the monokaryon or dikaryon (hybrid) P5 transformants. 
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Test crossing and fruiting 

To verify the biologically active transformed genes, transfor- 
mants were crossed with monokaryotic tester strains harboring 4 
different mating types (Table 1). Crosses between monokaryotic 
cultures were performed by placing actively growing mycelial 
plugs approximately 1 cm apart in the center of a Petri dish of 
MCM [12]. The mycelia were screened for clamp connections 
under 40 x magnification at the growing boundary lines of either 
side of the interacting strains. Fruiting tests for all crosses were 
performed as described by Ryu et al. [12]. 

Results 

Isolation and sequence analysis of the mating 
type-specific region 

The B3 mating-type locus was isolated using PGR amplification 
with the fi3-specific SCAR primer sets from a genomic library 
(more than 8,850 fosmid clones) constructed from KNR2312. 
DNA inserts of P. eryngii in cPE932-22 and cPE1462-8 were 
subcloned and sequenced using the shotgun method. Only the 
subclones containing 2-4 kb of DNA were sequenced in the 
forward and reverse directions using universal sequencing primers. 
The sequence data were edited, aligned, and merged using 
software. The sequence sizes in cPE932-22 and cPE 1462-8 were 
40.4 kb and 32.7 kb, respectively, and the total size of the 
sequences in the two clones was approximately 60 kb, excluding 
overlapping sequences. The SCAR marker [12] is located at the 
center of the locus (Figure 1), where the sequences from the two 
clones overlapped; thus, the two fosmid clones likely covered the 
complete B3 locus. 

Identification of the pheromone and pheromone 
receptor genes on the B3 locus 

ORFs and their putative functions were identified in the entire 
60-kb B3 locus using BLASTX and FgeneSH (Softberry). A 
physical map of the B3 locus is shown in Figure 1. The four 
pheromone and four receptor genes in the B3 locus of KNR2312 
spanned less than 12 kb. To verify the existence of subloci in 
KNR2312, flanking sequences were subjected to TBLASTX and 
contained adcUtional open reading frames, but none of the 
encoded proteins showed sequence similarity to known phero- 
mone or receptor genes (Table 2). However, we did not exclude th 
possibility that there are B3 subloci in KNR232 even though there 
was no non-parental B locus in the 98 meiotic monokaryons from 
KNR2312 (12). 



The CLA4, AKORl, and AKOR2 genes, which are located 
next to or near the B mating locus in Pleu7vtus djamor [14], were 
not identified at —26.3 to +21.7 kb relative to the B3 region 
(Table 2). Instead, cytochrome P450 {CYP) and reverse transcrip- 
tase-RNase H-integrase (pol) was identified upstream of the B3 
region. 

All receptors contain a 7 -transmembrane (TM) motif with 
varying amino acid lengths: 328 for PESTE3.3.1, 420 for 
PESTE3.3.2, 558 for PESTE3.3.3, and 470 amino acids for 
PESTE3.3.4. Alignment analysis revealed that the identities 
between PESTE3.3.1 and PESTE3.3.2 were 64.8% (full length) 
and 83.4% (compared with the truncated PESTE3.3.2 [326 amino 
acids], which did not include the non-corresponding C-terminal 
sequence of PESTE3.3.1). Comparison analysis showed that the 
C-terminal cytoplasmic domains varied in length across the genus, 
species, and even within a locus. PESTE3.3.1 has only 34 amino 
acid residues in its C-terminal cytoplasmic domain (Table S3). The 
C-terminal cytoplasmic domain of 5. commune was relatively long 
(326-344 amino acids), whereas those of U. rmydis and Lentinus 
edodes were 54 and 26-35 amino acids long, respectively. 

Using maximum likelihood algorithms, the molecular phyloge- 
netic analysis of receptors from P. eryngii and several mushroom 
fimgi, which are available in GenBank (http://www.ncbi.nlm.nih. 
gov/ genbank/), showed that two distinct clades were grouped with 
basal STE-like proteins such as U. rmydis PRA2 (Figure 2). Clade 
2 was roughly sub-divided into 3 groups, with not aU having strong 
bootstrap support (>80). PESTE3.3.1 and PESTE3.3.2 are part of 
clade 2, which were similar to P. ostreatus D330 STE3-1 and P. 
eryngii D625 STE3, whereas PESTE3.3.3 belongs to clade 1 with 
L. hicolor LbSTE3.2 and C. cinerea RCB2.42. PESTE3.3.4 
belonged with the receptors of clade 1 along with the orthologous 
C. cinerea RCB2.43 and S. commune BAR8. Sub-clusters 
containing P. eiyngii receptors possessed a known mating-type 
gene, indicating that the receptor genes from P. eryngii likely 
evolved from the common ancestor of the known mating genes. In 
particular, L. edodes receptors were found in a basal position in 
every subgroup of clade, suggesting that the sequences of L. edodes 
receptors are most likely closer to the ancestral gene than the other 
receptor pheromone sequences in clade 2. 

Fungal pheromone genes are sufficiently different to not 
hybridize with other genes, but they share consensus sequences 
in the C-terminal end of the protein precursors, the CaaX motif in 
C. cinerea and the CXXX motif (CVCH, CVRG, CWA) in S. 
commune [9]. CaaX is conserved for post-translational modifica- 
tion in most pheromones; therefore, it can be used to identify the 
clade of the pheromone. We translated all possible peptides and 
evaluated the CaaX sequence at the C-terminal end. In the P. 



-fj- 



CPE932-22 



CPE1462-8 



-fj- 



PESTE3. 3. 1 PESTE3. 3. 2 PESTE3. 3. 3 



PESTE3. 3. 4 



I^phbS.l PEpU)3.2 PEphb3.3 



PEphb3. 4 



Ikb 



Figure 1 . Physical map of tfie B3 locus. The gDNA sequence of the 63 locus is covered by two overlapping fosmids, cPE932-22 and cPEI 462-8. 
The genes identified in the sequences are shown as a horizontal wide arrow. The arrows indicate the direction of transcription. The gray bar on the 
map represents the sequence amplified by the 63 locus-specific primers (13-22ioo K F). PESTE3.3.1 stands for Pleurotus eryngii STE3 at the 63 locus and 
2" (serial number) gene, PEphb3.1 stands for Pleurotus eryngii pheromone at the 63 locus and V'' (serial number) gene. 
doi:1 0.1 371/journai.pone.01 04693.g001 
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Table 2. Gene homologs identified in tlie P. eryngii B3 mating-type locus and tlieir flanl<ing regions. 





Gene 


£ value^ 


Position 


Homolog 


Possible function 


rluF 


4.3 


245-2109 


ribosomal large subunit pseudouridine synthase F 
[Flavobacteria bacterium BBFL7] (ZP_01201616.1) 


responsible for the synthesis of pseudouridine 
from uracil-2604 in 23S ribosomal RNA 


Pol 


28-" 


3536-4072 


DNA/RNA polymerase, partial [Trametes versicolor 
FP-101664 SSI] (EIW61623.1) 


DNA repair 


Pol 


38-'' 


5603-6601 


reverse transcriptase-RNase H-integrase 
[Triciioloma matsutal<e] (BAA78625.1) 


DNA integration and RNA-dependent DNA 
replication 




1.7 


9893-10651 


conserved hypothetical protein [Trichophyton 
verrucosum HKI 0517] (XP_00301 8684.1) 


unknown 


CYP 




13695-15784 


cytochrome P450 [Phanerochaete chrysosporium] 
{BAL05081.1) 


catalysis of the oxidation of organic substances 


omt 




16301-23699 


0-methyltransferase [Punctularia strigosozonata 
HHB-11173 SS5] (EIN06883.1) 


catalysis of the transfer of a methyl group to the 
oxygen atom of an acceptor molecule 


rnt 


0.03 


24383-25349 


ribonuclease T [Bermanelta marisrubri] 
{ZP_01305716.1) 


catalysis of the endonucleolytic cleavage into 
nucleoside 3'-phosphates and 3'- 
phosphooligonucleotides ending in Gp with 2',3'- 
cyclic phosphate intermediates. 


PEphb3. 1 


6.6 


26325-26510 


fungal mating-type pheromone [Coprinopsis 
cinerea okayama7#130] (XP_002910434.1} 


mating type-specific peptide pheromone 


PEphb3.2 


1.7 


27056-27232 


B mating type pheromone precursor [Ceriporiopsis 
subvermispora B] (EMD35227.1) 


mating type-specific peptide pheromone 


PESTE3.3.1 


6e-"^ 


28186-29384 


putative pheromone receptor 
[Flammulina velutipes] (AGG68705.1} 


mating type-specific G-protein-coupled receptor 


PEphb3.3 


9.7 


29629-29877 


fungal mating-type pheromone [Coprinopsis 
cinerea okayama7#130] (XP_001 83497.1) 


mating type-specific peptide pheromone 


PESTE3.3.2 




30416-31911 


putative pheromone receptor 
[Flammulina velutipes] (AGG68705.1) 


mating type-specific G-protein-coupled receptor 


PESTE3.3.3 




32405-34232 


hypothetical fungal pheromone GPCR, STE3-type 
[Postia placenta Mad-698-R] (XP_002471247.1) 


mating type-specific G-protein-coupled receptor 


PEphb3.4 


0.46 


35693-35869 


lipopeptide mating pheromone precursor bap3-1 
[Schizophyllum commune H4-8] (XP_003028251.1) 


mating type-specific peptide pheromone 


PESTE3.3.4 




36579-38224 


pheromone receptor Rcb2 B43 [Coprinopsis 
cinerea okayama7#130] (XP_002910431.1) 


mating type-specific G-protein-coupled receptor 




18-'" 


38948-41569 


hypothetical protein SERLA73DRAFT_78193 [Serpula 
lacrymans var. lacrymans S7.3] (EGN93817.1) 


unknown 


ubIA 


3e-'° 


43157-43972 


UbiA prenyltransferase [Saprospira grandis str. 
Lewin](YP_005323559.1) 


integral to membrane 




9e-=° 


44187-45557 


predicted protein [Laccaria bicolor S238N-H82] 
(XP_001 881 105.1) 


unknown 


sgnH 


2 e-™ 


46375-48010 


SGNH hydrolase [Trametes versicolor FP-101664 
SS1](EIW56391.1) 


lipid metabolic process 




0 


48339-51782 


hypothetical protein SERLA73DRAFT_1 72732 [Serpula 
lacrymans var. lacrymans S7.3] (EGN92622.1) 


unknown 




0 


53063-55751 


hypothetical protein SERLA73DRAFT_1 90801 [Serpula 
lacrymans var. lacrymans S7.3] {EGN92621.1) 


unknown 


Pol 


0 


57385-59974 


reverse transcriptase-RNase H-integrase 
[Tricholoma matsutake] {BAA78625.1) 


DNA integration and RNA-dependent DNA 
replication 



^The expected value (E) indicates the number of hits one can "expect" to see by chance in GenBank similarity searches. Under "Homolog", the species form with the 
lowest f-value hit or the value closest to zero is given in brackets, followed by the GenBank accession number in parentheses. 
doi:1 0.1 371 /journal.pone.Ol 04693.t002 



eryngii B3 locus, four pheromone precursors were identified using 
DNAMAN (Lynnon Biosoft, Cambridge, UK) and BLASTP. A 
BLASTP search sliowed that each pheromone precursor has a 
similar precursor, including C. cinerea phb3.2 in B43 
PCP_002910434.1) for PEphbl.l, the Ceriporiopsis subvermispora 
B mating-type pheromone precursor (EMD35227.1) for 
PEphb3.2, C. cinerea fongal mating-type pheromone 
{XP_00 1834397.1) for PEphb3.3, and phb3 in the B5 locus in 
C. cinerea (AAQ96360.1) for PEphb3.4. The 4 pheromone 



precursors contain 57-61 amino acids and have a cysteine residue 
at the 4th site from the C-terminus. During post-translational 
modification, 3 amino acid residues at the C-terminus are 
removed, and a 4th cysteine residue is farnesylated. Mature 
pheromones were 9-15 amino acids in length [8,21] in the B3 
locus; the mature pheromone size was assumed to be 12 to 13 
amino acids based on the consensus sequences for the excision site 
(Figwe 3). Additional consensus sequences were found at the ER 
or TR at the 9-1 5th sites from the C-terminus and at a proline at 
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Figure 2. Phylogenetic relationships of the fungal pheromone receptor protein sequences. CC: C. cinerea; SC. S. commune. rcb1, rcb2, and 
rcb3 from C. cinerea; Ba. and 6(3 from S. commune; PDSTE3.3 from P. djamor; PESTEB-like from P. eryngii; and LBSTE3-like from L. bicolor. The GenBank 
accession numbers of the receptors analyzed here are shown in Table 2 and were previously reported [14], except the following; L edodes SUP2 B2 
rcb2: AER51012.1; V. volvacea ste3.1: JX982139.1; V. volvacea ste3.2: JX982140.1; V. volvacea ste3.3: JX982141.1; V. volvacea ste3.4: JX982142.1; F. 
velutipes W23 ste3.2: KC208605.1; F. velutipes W23 ste3.4: KC20861 1.1; F. velutipes KACC42780 ste3.1: HQ630590.1; and F velutipes KACC42780 ste3.2: 
HQ630591.1. 

doi:10.1371/journal.pone.0104693.g002 



the —4 and aspartate/ asparagine at —2 relative to the cutting site. 
These conserved sites were recognized by processing molecules 
[9]. Sequences are typically more diverse across a genus and 
species than within a species. The cDNA sequences of the 
pheromone and receptor genes were deposited in GenBank with 
the primers used for their amplification. 

Gene expression of the pheromone receptor and 
pheromone genes in monokaryons and dikaryons 

Real-time qRT-PCR was performed with KNR2312 (dikaryon) 
and P6 (monokaryon) to quantify the expression of the 4 
pheromone precursor genes and the 4 receptor genes. Gene 
expression was normalized to 18S rRNA expression and calibrated 
to the assigned value of 1 by the standard curve method (Qiagen). 



The transcription level of the receptor genes was significandy 
higher (95% confidence interval) in the monokaryon (P6) (3.8-12.2 
fold) than in the dikaryon (KNR2312), except for PESTE3.3.2 
(Figure 4a, Table S4). The expression patterns of the pheromone 
precursor genes were similar to those of the receptor genes, 
ranging from 1.5-31.4 fold (Figure 4b, Table S4). In particular, 
PESTE3.3.4 and PEphb3.2 showed the most significandy 
different expression pattern between moiiokaryotic and dikaryotic 
mycelia. 

Function of pheromones and pheromone receptors in 
the fruiting process 

The transgenic monokaryon P5 was developed using two 
pheromone genes and one pheromone receptor gene by 



Name size 

P. eryinffii PEphbS. 1 HDTFDLLETFILSDSESESSSCMVITDIPSSAEODADVFDLPIH HERDVPGGTFRGWC I I S 60 

P. eryrig-ii PEphbS. 2 MDTFYTITAITDALEQCIAPVLSSIPIQNMSDDDPPLPVPVDE ERYGGGVTSSWC TIS 57 

P. erroffii PEphbS. 3 MDTFMSLDALEQLISTILPPSISILIPIESHACCIDDDELPIPIHE ERYGGGGVTSSWC IIS 61 

P. eryngii mpH>3.i MDTFISLDEFFVEGSASDVDIMSSSDSEGEDSDCEAPMPIDF ERRGPGGAAOFWC IIS 57 

C. cinerea okayamaT MDSFSSLDAFWDEELLOLAIDIPOPIPGODOOPPINE ERPGAGTMGAFC IIH 52 

C. siAvermsporai MDAFFTLAAPVPSADETEIQVPVDE EQFGSSTTSQHC IIA 39 

Figure 3. Alignment of the pheromone precursor amino acid sequences. The processed pheromone sequences are underlined, and the 
putative C-terminal doublets (ER) are shown. The GenBank accession number of the C. cinerea okayama 7 pheromone is XP_001 834397.1, and that of 
Ceriporiopsis subvermispora B is EMD35227.1. 
doi:1 0.1 371/journal.pone.01 04693.g003 
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Figure 4. Relative expression of P. eryngii KNRIS^ 2 receptors 
(a) and pKieromones (b) in the monoicaryon (P6) and diluryon 
(KNR2312), as determined by real-time qRT-PCR. Gene expres- 
sion was normalized to 18S rRNA expression and calibrated to the value 
for the monokaryon (P6), which was assigned a value of 1, by the 
standard curve method (Bio-Rad). All assays were performed in 
triplicate. The error bars show the standard deviations for triplicate 
samples. 

doi:10.1371/journal.pone.0104693.g004 



Agrohacterium-mediMed transformation. Twelve transformants of 
P5phb31, 16 of P5phb33 and 3 of P5STE331 were screened using 
hygromycin selection medium and were confirmed by PGR with 
A^/j-specific primers. PGR with the SSR45 primer set produced 
one single band identical to that of P5, indicating the monokary- 
otic nature of the transformants (Figure 5a). Glamp connections 
were observed in dikaryon crosses between three transformants 
and KNR2312-2 {A4B4) (Figure 5b and 5c (3)-(5)), whereas 
crosses between transformants and monokaryons [A3B3 or A3B4) 
could not develop clamp connections due to incompatible 
combinations (Figure 5b). The control combination of P5 and 
the tester strain {A4B3) possessed clamp connections (Figure 5c 
(1)), whereas the transformant containing the pBGgHg vector was 
unable to form clamp connections when crossed with the A4B4 
tester strain (Figure 5c (2)) due to incompatibility at the A locus. 
The positive control developed the clamp connections within a 
week, whereas the transformed monokaryons required 30-40 days 
after the two mycelia contacted. Moreover, the density of the 
clamp connections was 2.5 times higher in transformed mono- 
karyons than in the wild-type strain (Figure 5c). Nonetheless, on 
the basis of our results, the transformed pheromone and receptor 
genes appear to play a role in the clamp connection process. 

AH combined dikaryons, including the positive and negative 
control strains, were fruited in sawdust medium. Only compatible 
combinations showed fruiting bodies similar to those of the 
positive strain (Figure 5c). The dikaryons P5phb31 xKNR2312-2 
and P5STE331 XKNR2312-2 showed full fruiting bodies, whereas 
those of P5phb33xKNR2312-2 did not fully differentiate and had 
primordia without caps. All fruiting bodies were relatively small 



and delayed (more than 25 days) compared with KNR2312 (18 
days, data not shown). 

Discussion 

Structure of the B3 locus 

Mating-type genes show extraordinarily diverse sequences; thus, 
using sequence-based methods to isolate whole mating-type genes 
is difficult. Although whole genome sequencing has been 
widespread and most of the genes of interest are available in 
online databases, the pheromone and receptor genes have distinct 
sequences in each strain. Alternatively, genomic subtraction and 
positional cloning using mitochondria intermediate peptidase 
(M/P), CLA4, putative multidrug transporter (mfsl), and PABl 
have been employed to isolate the mating region [5,8,14,22]. 
Additionally, marker genes are sometimes located too far from the 
B mating region to isolate the entire sequence of interest [14]. 
Recently, we identified mating-type alleles in the A and B loci and 
developed SGAR primers specific to the B3 locus. Moreover, the 
SGAR marker, shown to be located in the middle of the B3 locus, 
was amplified in AxB3 monokaryons from in 100 monokaryons 
[12]. We isolated the B3 mating region and identified a variety of 
genes, including mating-type genes, by sequence analysis using 
Fgenesh, BLASTX, and manual curation. The B3 locus of P. 
eryngii possesses one more receptor (four receptors) than the 
typical B locus in C. cinerea (three receptors). In S. commune, the 
B locus is divided into two subloci, 5a and 5/? [19,20], whereas in 
C. cinerea, only a single B locus is present. More than three 
pheromone genes at the B locus have been reported in C. cinerea, 
Flammulina velutipes, and L. hicolor. In most of these cases, unlike 
the B3 locus of P. eryngii, extra receptors (other than the three 
functional receptors) that did not co-segregate with the mating 
type were located outside the locus and were considered to be 
impaired remnants of duplications [4,15]. According to the 
phylogenetic analysis, PESTE3.3.1 and PESTE3.3.2 appear to 
be very closely related (Figure 2), and the sequence identity (64.8% 
between intact sequences and 83.4% between truncated and intact 
sequences) suggests that these two proteins were most likely 
duplicated later than any other receptors in the B3 locus. The 
expression of the receptor and pheromone genes was supported by 
real-time qRT-PGR data. However, due to the short length of the 
G-termiiial cytoplasmic residues, PESTE3.3.1 seems likely to be a 
receptor-like gene. Whether such gene types play a role in the 
interaction between the receptor and downstream molecules 
remains unclear. Although P5-STE3.3. 1 xKNR2312-2 showed 
sparse clamp connections and late fruiting compared with the 
positive control (P5xKNR2312-6), it produced similar fruiting 
bodies to that of the positive control, indicating that a receptor 
with short G-terminal cytoplasmic residues could play a role in the 
fruiting process. In the case of S. commune, the fi-dependent 
mating reaction was limited and produced an abnormal fruiting 
body after the shortening of the intracellular G-terminus [23]. In 
contrast, the truncation of the G-terminal intracellular region of 
the S. commune receptors has not been observed to affect signal 
transduction [24-26]. Furthermore, this short cytoplasmic region 
was found in not only P. eryngii but also other mushrooms such as 
U. tnaydis and L. edodes (Table S3). In addition to the C-terminal 
tail, multiple sites including the intracellular loop domains in 7- 
TM are thought to be involved in mating-type specificity 
[24,25,27]. Further analysis is required to determine how the 
length of the cytoplasmic domain of the receptor functions in the 
mating process. Novel pheromone-receptor interactions presum- 
ably evolved through gene duplications of the prototype STE3-like 
gene and the pheromone gene and subsequent mutation of the 
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Figure 5. PCR band pattern for the confirmation of transformation, crossing table between the wild-type and transgenic 
monokaryons, and fruiting body (clamp connection) formation, a; Specific DNA fragments obtained in test-crossed mycelia using hph (Table 
SI). P: positive marl<er (pBGgHg), KNR2312: a dil<aryon consisting of P5 and P6, P6: a monol<aryon, P5: a monol<aryon, pBGgHg: a transformant with 
the empty vector (pBGgHg), P5phb31: a P5 transformant with PEphbS.l, P5phb33: a P5 transformant with PEphb3.3, P5STE331: a P5 transformant with 
PESTE3.3.1, b: Wild-type tester strains (KNR2312-4: A3B3, KNR2312-6: A3B4, KNR2312-29: A4B3, and KNR2312-2: A4B4) were used for crosses with the 
transgenic monokaryons with the empty vector or with the pheromone or receptor genes, c: Clamp connections and fruiting bodies were observed 
in the positive control (wild type) (1) and dikaryons mated with transgenic monokaryons and the tester strain (3 and 5), but none were observed in 
the incompatible combination (2). In the case of 4, only primordia were shown. Labels above the panels indicate the cross combinations. 
doi:1 0.1 371 /journal.pone.01 04693.g005 



duplicated genes (for review see [28]). Because single- or double- 
amino acid substitutions in receptors or pheromones can 
significantly alter interaction specificity [29], PESTE3.3.1 could 
be a candidate for new mating specificity. 

Gene expression of the pheromone receptor and 
pheromone genes in monokaryons and dikaryons 

Few studies have been conducted to analyze the expression of 
mating-type genes in mushroom fungi. Monokaryon and 
dikaryon fungal states show different expression levels of the A 
and B loci genes. Richardson et al. observed constitutive 
expression of the A locus genes to a similar extent in both 
monokaryons and dikaryons of C. cinereus [30]. In another 
study [23], a different expression pattern was observed for the 
receptor and pheromone genes of S. commune; at 12 h (mid- 
point of massive dikaryosis) after mating, both genes were highly 
expressed compared with the initial stage; at 72 h (end of 
massive dikaryosis), the pheromone genes were expressed at 2- 
fold higher levels than in the monokaryotic state, whereas the 
receptor genes were expressed at similar levels in the mono- 



karyotic and dikaryotic states. The truncated pheromone was 
expressed at a lower level in the dikaryon. In the present study, 
comparative expression of the receptors and pheromone 
precursor genes revealed that most pheromone precursor genes 
were highly expressed in the monokaryons compared with the 
dikaryons (Figure 4), indicating that the fungal monokaryon 
status requires additional mate-specific pheromone and receptor 
genes for a monokaryon-monokaryon interaction to generate a 
dikaryon. The roles of pheromone perception in Saccharomyces 
cerevisiae have been described previously in mate attraction and 
the formation of shmoo cells that grow towards their mating 
partner [31]. Similarly, in Neurospora crassa, pheromones and 
receptors were involved in the chemotropic polarized growth of 
female-specific hyphae (trichogynes) toward male cells of the 
compatible mating type [32,33]. A dikaryon does not require 
this attraction towards a compatible gamete. Before the mating 
process, sufficient amounts of the pheromones and receptors are 
required because the opposite mating type can only be met 
randomly. 
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Function of pheromones and pheromone receptors in 

the fruiting process 

Fungal transformation, especially in edible mushrooms, is not 
easily achieved, although some methods have been developed 
[34,35]. We adapted Agrobacterium-mtAmteA transformation, 
thus succeeding in the isolation of several transformants. The 
dikaryons from crosses between transformants {A3B4::PEphb3.1 , 
A3B4::PEphh3.3, A3B4::PESTE3.3.1) and a monokaryotic 
tester strain {A4B4) demonstrated clamp connections, suggesting 
that the transgenic genes most likely produced compatible 
combinations with the tester strain. Although clamp connections 
were formed more slowly and sparsely in the transgenic dikaryons 
than in the positive control (Figure 5c), the transgenic dikaryons 
showed full fruiting bodies, except P5phb33 xKNR2312-2, which 
indicated the requirement for more components of the normal 
process. This result is consistent with a previous report [23], in 
which malformed fruiting bodies and spores were observed in the 
transgenic dikaryons. The small and delayed fruiting bodies were 
likely caused by selfing depression. 

Supporting Information 

Table SI Primers for the isolation of the B3 locus and 
the confirmation of transformation and expression. List 
of specific primer sets for the isolation of the B3 locus from the 
gDNA library, the confirmation of transgenic mycelia and the 
amplification of pheromone and receptor gDNAs with their 
promoter and terminator sequences. 
(DOCX) 
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